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Abstract 

It is shown, that highly accurate estimation of muon shower content can be performed on the basis 
of knowledge of only vertical depth of shower maximum A^ ax and total signal in ground detector. The 
estimate is almost independent on primary energy and particle type and on zenith angle. The study is 
performed for 21500 showers, generated with CORSIKA 6.204 from spectrum E~ x in the energy range 
log 10(E) [eV]=18.5-20 and uniformly in cos 2 6 in zenith angle interval 6 = 0° -65° for QGSJET II/Fluka 
interaction models. 

Muon shower content is the key parameter for studies of primary mass composition and test of hadronic 
interaction models. Using universality property of extensive air showers (EAS) development [1-6] we propose 
a simple and precise method to determine muon shower content from vertical depth of shower maximum 
and total signal (signal in water tanks or in scintillator detector). The study is performed for 21500 showers, 
generated with CORSIKA 6.204 [7] from spectrum E~ x in the energy range logl0(£') [eV]=18.5-20 (with 
different statistics in 3 energy bins 18.5-19.0, 19.0-19.5 and 19.5-20.0) and uniformly in cos 2 8 in zenith 
angle interval 9 = 0° - 65° for QGSJET II [8-10]/Fluka [11, 12] interaction models. Electromagnetic (EM) 
component thinning was set to 10 -6 , the observation level was at 870 g/cm 2 . All longitudinal showers 
characteristics and charged particles density were taken directly from CORSIKA output files. The expected 
signal in Auger-like tanks was calculated according to the procedure described in [13, 14] with the use of the 
same GEANT 4 lookup tables as in [14]. 

From different aspects of universality of shower development we will be interested only in the dependence 
of electromagnetic and muon signals on the distance of shower maximum to the ground. We will begin 
our consideration from Auger-like experimental setup and consider signal in water Cherenkov tanks at 1000 
meters. In this case one deals with the situation, where muon contribution to the detector is much larger 
compared to the EM one due to higher energy losses of muons in tanks. The common way to express 
universality of electromagnetic signal is to plot it against slant distance to the ground (Fig.fl] and also 
Figure 1 in [4]), showing its quasi- independence on primary particle type. The muon signals functional 
dependence on slant distance to the ground DG is also very similar for both proton and iron, but there 
is a shift in the normalization (Fig. [I]). Since iron showers reach X max earlier than showers from protons, 
comparison of set of showers from p and Fe at equal DG means the comparison between showers with 
different zenith angles, but at the same development stage. Passing to comparison of shower characteristics 
dependence on vertical distance to ground DG V reveals a very interesting property (see Fig. [2]): in this 
case the similarity of functional dependence of muon and EM signals on DG V between p and Fe primaries 
is preserved, but now also EM signal normalizations are different. This happens because one confronts 
showers, which have the same vertical distance from A^ ax to the ground, but once again proton showers 
are more inclined than iron ones and their EM component attenuates more while reaching the ground from 
I^ M . The ratio S'em/'S'em turns ou t to be almost equal to the S^/S? one and this allows to state the 
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new shower universality property: the ratio of the muon signal to the EM one S^/Sem is the same for all 
showers, reaching the maximum at the same vertical depth X* m , independently on the primary particle type, 
primary energy and incident zenith angle (at the least for the energy and angular ranges considered here). 
This property is illustrated in Fig. [31 where it is shown the dependence of S m /Sem on X^ ax for p, O and Fe 
primaries in four different energy bins. The functional dependence between X^ mK an d S^/Sem turns out to 
be very simple and quasi-universal for all energies and primaries. The function in the form 

= M S n/ Sem + a) h (1) 

fits well the data and the fit parameters are quite stable across entire energy range. Having in hand the 
functional dependence of A^ ax on (S^/Sem) and using S to t — Sem + one easily gets the equation, which 
allows to obtain muon signal from shower vertical depth and total signal in tanks: 

gfit _ gtot 

" l + l/{{X^JA) 1/b -a) 

We calculated the difference between the Monte-Carlo (MC) simulated muon signal S^f and the muon 
signal, obtained from the fit Sftf, an example distributions of this value are shown in Fig. [4j In Table Q] we 
give mean and RMS values of such distributions for various energy bins, obtained with the unique set of fit 
parameters A = 538, b = —0.25 and a = —0.22. It is seen, that the estimates of muon signals are unbiased 
with less than 1% deviation of mean reconstructed muon signal from the MC one for all 3 primaries and 
the RMS values are small: 8% for protons and around 5% for oxygen and iron. Certainly, application of 
specific coefficients for every energy bin or narrowing of zenith angle interval, or using of more sophisticated 
fit functions can even slightly improve the performance of the method, which anyhow is good in its simple 
and universal form. The described universal dependence of S^/Sem on X^ ax holds true in the wide interval 
of distances and in Fig. [5] we show examples for the distances 200 and 1500 meters from the core, though for 
distances closer to the core the function in the form (fT]) does not describe accurately the data in the entire 
angular range 0° — 65° and it is needed or to split it in two parts or to apply more complex parametrization. 

The same universality principle holds true in the case of detectors using scintillators and effectively 
measuring density of charged particles [15]. We have performed the reconstruction of muon densities using 
the dependence of the ratio (muon density D^)/ (electron density D e ) on X^ &K in the form same to (p} (see 
Tableland Fig. [6]). It is seen, that when muons and electrons equally contribute to the detector signal 
the shower fluctuations play more important role and the accuracy of parametrization is only within 15%, 
though the estimate is still unbiased. 

Hence, the new universality property allows to obtain accurate estimates of the muon signal, which are 
almost independent on the primary particle type, primary energy and zenith angle for various types of 
ground detectors. Taking in consideration that the shower universality property was established for different 
interaction models [1-6], we expect that the proposed approach to muon content derivation is not specific 
only to QGSJET II/Fluka case. 

The discovered simple shower universality property in respect to (muon signal/EM signal) ratio, giving 
access to the muon shower content, can open new possibilities as in solution of the global problems, such 
as derivation of primary mass composition and understanding of hadronic interactions properties, so in the 
number of more particular tasks (e.g. estimations of primary energy on the basis of pure electromagnetic 
signal, primary particle type independent corrections to the missing energy in experiments using fluorescent 
light etc.). 
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Figure 1: EM and muon signals from proton (red squares) and iron (blue crosses) in water Cherenkov tanks 
at 1000 m vs slant distance from shower maximum to the ground DG in log 10(1?) [eV]=18. 9-19.0 energy bin 
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Figure 2: EM and muon signals from proton (red squares) and iron (blue crosses) in water Cherenkov tanks 
at 1000 m vs vertical distance from shower maximum to the ground DG V in log 10(-E)[eV]=18. 9-19.0 energy 
bin 
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Figure 3: Ratio of signals in water Cherenkov tanks S^/Sem at 1000 m vs vertical depth of shower maximum 
X^ ax in four energy bins. Protons — red squares, oxygen — brown diamonds, iron — blue crosses 
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Figure 4: Distributions of relative difference between MC simulated muon signals in Cherenkov water tanks 
S^ G and muon signals derived from the fit S** at 1000 m. Protons — red line, oxygen — brown line, 
iron — blue line. 
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Figure 5: Ratio of signals in water Cherenkov tanks S^/Sbm at 200 m and f500 m vs vertical depth of 
shower maximum X^ ax in log 10(£')[cV]=18. 9-19.0 energy bin. Protons — red squares, oxygen — brown 
diamonds, iron — blue crosses 
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Table 1: Means and RMS of distributions of relative difference between MC simulated muon signals in 
Cherenkov water tanks S^f c and muon signals derived from the fit at 1000 m (see also Fig. 0]) (S^f G — 
S , ^ t )/S'^ IC , %, calculated with the unique set of parameters for all energy bins: A = 538, b — —0.25, 
a = -0.22 
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8.3 


-0.3 
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-0.2 
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0.4 
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0.2 


5.1 


-0.1 
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Figure 6: Top: ratio of muon density to the electron one at 1000 m vs vertical depth of shower maximum 
^Qnax f° r two energy bins; bottom: distributions of relative difference between MC simulated muon density 
-D^ 1 and muon density derived from the fit at 1000 m. The data are given for two energy bins, 
protons — red, oxygen — brown, iron — blue. 
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Table 2: Means and RMS of distributions of relative difference between MC simulated muon density Djf 
and muon density derived from the fit Df at 1000 m (see also Fig. EJ) {Df c - D^)/D^ c , %, calculated 
with the unique set of parameters for all energy bins A = 475, b = —0.28, a = —0.09 
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